The McVittie solution with a negative cosmological constant 
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Whereas current cosmological observations suggest that the universe is dominated by a positive 
cosmological constant (A > 0), the AdS/CFT correspondence tells us that the case A < is still 
worthy of consideration. In this paper we study the McVittie solution with A < 0. Following a 
related study, the solution is understood here by way of a systematic construction of conformal 
diagrams based on detailed numerical integrations of the null geodesic equations. As in the pure 
Robertson - Walker case, we find that A < ensures collapse to a Big Crunch, a feature which 
completely dominates the global structure. 
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I. INTRODUCTION 



written in the form (e.g. Q) 



Recently [l| , a detailed study of the McVittie solution 
2] was carried out for a non-negative cosmological con- 
stant (A > 0). The McVittie solution has been known 
for many years, but it continues to attract interest Q. 
Even though it is now widely believed that the universe 
is dominated by a positive cosmological constant, the re- 
markable AdS / CFT correspondence [4j presents a strong 
argument that the case A < should also be examined. 
Following [lj we systematically construct a global view 
of the McVittie solution with A < based on numeri- 
cal integrations of the null geodesies. What results is a 
situation very distinct from the case A > 0: the global 
structure is completely dominated by a collapse to a Big 
Crunch, just as in the pure Robertson - Walker case. 



II. THE SOLUTION 

A. Overview 

For a perfect fluid with energy density p and isotropic 
pressure p the strong energy condition [5[ is given by 



p + 3p > 0. 



(1) 



For a Robertson - Walker background with scale factor 
a(t), for (TTJ, Einstein's equations with A < give 



ds 2 



1 — m/2u 
1 + m/2u 



dt 2 + a 2 (1 + m/2u) 4 (dr 2 + r 2 dVL 2 2 ) 

(3) 

where u = ra, m is a positive constant and dfl 2 is the 
metric of a unit 2-sphere. If the McVittie solution ^ 
asymptotes to a reasonable Robertson - Walker back- 
ground, then, as in the pure Robertson - Walker case, 
A < gives rise to a Big Crunch which dominates the 
global structure. As previously [l[ , we use the coordinate 
transformation 



to obtain 



R(t,r) = u(\ + m/2uf 



2H(t)R 



(4) 



dR 2 



ds 2 = -f{t,R)dt 2 v > dtdR+ z — —+R 2 dfl 



where 



y/1 - 2m/R 1 - 2m/ R 



f=l- 2m/R-H 2 R 2 



(5) 



(6) 



and H is the Hubble function a /a. From ([5]) it follows 
that tangents to surfaces of constant finite t are spacelike 
for R > 2m (and so for finite t we set the future orien- 
tation dt/dX > for affinc A increasing to the future) 
and tangents to surfaces of constant R are spacelike for 
/ < 0, null for / = and timelike for / > 0. We note 
again that the effective gravitational mass Q associated 
with (|5|) is not m, but rather M, given by 



3- > -A > 
a 



(2) 



where ■ = d/dt, and so we necessarily have a Big Crunch 
The particular solution we are concerned with here 
is the simplest of the McVittie class, and this can be 
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M(t, R) = m + -H 2 R 3 . (7) 

B. The function H 

Quite unlike lj, we note that because of (0), 

* -> oo $. (8) 
Rather, we are now interested in models for which 

a(0) = a(tf) = 0. dfin) = (9) 



2 



where 



and 



< t < t f , 



(10) 



(11) 
(12) 



H(0<t< t a ) > 0, H(t ) = 0, H(to < t < t f ) < 
Note that from Q 

M(t ,R) = m, 
and from @ 

H < 0. (13) 

Whereas from the definition of it, limt-^o,*/ w = for 
all finite r, from the transformation ((4]) 



lim R 



if m = 
oo if m 7^ 



(14) 



As a result, neither t = nor t — tf are part of the 
spacetime j5) for ra ^ 0. We note that (|14p shows us 
that there is no continuous transition from m ^ to 
m = 0. 



C. Scalar Singularities 

As explained previously [l|, singularities, as revealed 
by scalars polynomial in the Riemann tensor, are re- 
flected here by the Ricci scalar TZ, 



K = 12H 2 + 



6H 



sjl - 2m/ R 



(15) 



since all other invariants, derived from (partial) deriva- 
tives of the metric tensor no higher than 2, add no 
new information in the cases under consideration. For 
< t < tf, since H ^ 0, there are singularities at 
R = 2m (u = m/2), which are spacelike. The appar- 
ent singularities at t = and at t = tj^\, over the range 
2m < R < co are, as explained above, not part of the 
spacetime. 



D. The locus / = 

As in the previous analysis fjj . the locus / = is im- 
portant for an understanding of the spacetime ([5])- How- 
ever, due to the nature of the function H studied here, 
this locus is quite distinct from the locus studied in [lj. 
In particular, the roots R± studied there do not exist for 



1 We assume that a ^ at t = and t - 



A < 0. First let us note that the locus f — includes 
R = 2m at t = to where H = 0. Moreover, since 



/ m 



H R) R = HHR 



(16) 



along the locus, R can have a vertical tangent in the R—t 
plane on < t < to and to < t < tf only at R — 3m. 
As R — > oo the locus becomes H 2 R 2 = 1 which requires 
H 0, that is, t -t t^- 



E. Null Geodesies - Qualitative 

The radial null geodesies of ([5|) satisfy 

yJl-2rn/R (HR ± y/l - 2m/ R) . (17) 



dR 
~dt 

We label the branch "+" "outgoing" and the branch 
"ingoing" . Clearly 



dR 
~dt 



= ± 1- 



2m 
~R 



(18) 



and so the ingoing geodesies already have dR/dt < at 
to whereas the outgoing geodesies have dR/dt > at to. 
Further, 



dR 
Hi 



= 



(19) 



at / = for ingoing null geodesies when t < t and for 
outgoing null geodesies when t > t . Since dR/dt > 
along both branches for / < and t < to and dR/dt < 
along both branches for / < and t > to it follows that 
the ingoing geodesies reach a maximum R at / = for 
t < to and the outgoing geodesies reach a maximum R at 
/ = for t > t - A special case is shown below in Figure 
□ and Figure H 



F. Energy conditions in general 

Let us start by rewriting the energy density and 
isotropic pressure in terms of H . From Einstein's equa- 
tions with A < we find 



8np = 3H 2 - A, 8tt P = -3H 2 - 



2H 



yji - 2m/ R 



A. (20) 



As a result, with (|13p . a general feature of these models 
is p > and p + p > and so the null and weak en- 
ergy conditions are always satisfied. The strong energy 
condition requires p + 3p > which, from (|2T)|) . gives 



H 



H 2 <±<0. 



yjl - 2m/ R 
The dominant energy condition requires 

- P < P < P- 



(21) 



(22) 
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It follows from ([13 ]) . (|20 |) and (122)) that whereas the left- 
hand inequality is generally satisfied, the right-hand in- 
equality requires 



H 



y/1 - 2m/ R 



< 3H 2 - A. 



(23) 



V. INTEGRATION OF THE NULL GEODESICS 



A. Integration in the R — T plane 



With (f25|) we can write the null geodesic equations (fi~7j) 
in the dimensionless form 



III. A SPECIFIC FORM FOR H 

We cannot proceed with further details without a spe- 
cific form for H. For notational convenience, and for a 
comparison with [l| , define 



(29) 

where Y = R/m and 5 is the parameter Horn. Numerical 
integrations of (|29[) are shown in Figures Q] and [2] 



Hi 



A 

-— T 

3' 



3Hot, 



(24) 



and take Hq > 0. For a we take the Robertson - Walker 
scale factor for spatially flat dust with A < 0. This gives 



H 



H sm(T) 
1 - cos(T) ' 



(25) 



Clearly H is periodic with period 2ir and so To = tt. 
Moreover, 



H = - 



<iH% 



1 - cos(T) ' 



(26) 



The dominant feature in our study is the development of 
a Big Crunch. We know that this will occur with A < 
as long as the strong energy condition holds. For more 
general equations of state, say p = np, the strong energy 
condition gives p(l + 3n) > 0. Since the minimum 8irp in 
our considerations is —A > 0, we would recover the same 
structure for k > — 1/3. In this sense our choice of k = 
is not critical. 




FIG. 1. Numerical integrations of (|29p and the locus / = 0. 
The ingoing geodesies (" — ") are shown dashed and reach 
their maximal value of R/m on the left branch of the locus 
/ = 0. The outgoing geodesies (" + ") are shown solid and 
reach their maximal value of R/m on the right branch of 
the locus f — 0. The future orientation is T increasing as 
shown. All geodesies begin and end at R = 2m which has two 
distinct parts, separated by the exceptional point T = tt, R = 
2m. This exceptional point is not part of the spacetime. An 
enlarged view near R = 2m is shown in Figure [5] 



IV. ENERGY CONDITIONS 



From (|2"Tj) we find that the strong energy condition is 
satisfied for 



y/l - 2m/ R 



> 2. 



(27) 



Since the left side of ([27)1 is at least 3, we see that the 
strong energy condition is always satisfied. For the dom- 
inant energy condition we rearrange ([23| to give 



R 8 
— > - 

m 3 



(28) 



and so the dominant energy condition fails sufficiently 
close to the singularities. 




FIG. 2. As in Figure [T] in the neighborhood of R — 2m. 
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B. Integration in the z — T plane 



VI. GLOBAL STRUCTURE OF THE 
SPACETIME 



As previously [H, we find it numerically convenient to 
compactify R and define 



2 
Y 



so that equation (|2"9|) takes the form 



(30) 



dz 
dT 



(1-z 2 ) ( sin(T) z{l-z 2 ) 



G 



1 - cos(T) 



± 



26 



(31) 



A. Construction of the conformal diagrarrQQ 

We represent the past boundary as a horizontal line in 
a Cartesian plane (y = 0, — f < x < 1), setting the right 
end of R = 2m at T = 0, and the left end at T = n. 
To represent the interval < T < tt, from x = — 1 to 1, 
we only require a one-to-one function between these two 
variables. Our choice for x(T) is given in the Appendix. 
This function was chosen purely for visual reasons by 
finding the spline curve fit of several points that were 
adjusted manually to make the resulting conformal dia- 
gram more visually appealing. Note that the choice of 
the shape of the curve (R = 2m, < T < tt) to be a hor- 
izontal straight line, as well as the function x(T) used is 
arbitrary and does not change the overall global structure 
presented in the conformal diagram. It is only required 
that the curve (R = 2m, < T < tt) be space-like, and 
the function x(T) be one-to-one. After this point, the 
procedure we used is identical to that in 



Numerical integrations of (1311) are shown in Figure [3] 
Any point in the spacetime (0 < T < 2ir, R > 2m) can 
be connected to the past boundary (R = 2m, < T < tt) 
by a unique null geodesic from each branch. We use this 
to construct the conformal diagram below. 



B. Null geodesies 

Under the procedure described above, Figure |3] is 
mapped into Figure |4j 





FIG. 3. Numerical integrations of (|31[) and the locus / = 0. 
The ingoing geodesies ( " — " ) are now shown solid and reach 
their maximal value of z on the left branch of the locus / = 0. 
The outgoing geodesies (" + ") are now shown dashed and 
reach their maximal value of z on the right branch of the 
locus / = 0. The future orientation is T increasing as above. 
All geodesies begin at z — for some value of < T < n and 
terminate again at z = for n < T < 2n. 



Conformal diagrams in the Robertson - Walker case can be found 
in HO]. 

Our understanding is that null afnnc distance is relevant iff the 
coordinates diverge, indicating the incompleteness of the coordi- 
nates for finite affine distances. There is no divergence of the co- 
ordinates in this paper. In particular, if the boundaries (R = 2m) 
are at finite null affine distance, this finiteness is irrelevant due 
to the fact that R = 2m is genuinely singular and no exten- 
sion is possible. If the boundaries (R = 2m) are at infinite null 
affine distance the conformal diagrams remain unchanged and 
complete. That is, the affine completeness of the null geodesies 
is, in this case, irrelevant, quite unlike the cases A > 
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FIG. 4. Conformal representation of the outgoing null 
geodesies (dashed) and ingoing null geodesies (solid). The 
locus / = is also shown. This locus is globally spacelike. 
Note that / > between the two branches of the locus. The 
boundaries are at R — 2m. The point at the lower right 
is R — > oo,0 < T < 2n. The point at the lower left is 
T = 7T, R = 2m. 



C. Surfaces of constant R and T 

Conformal representations of surface of constant R and 
constant T are shown in Figures [3] and [5] respectively. 




FIG. 5. The trajectories show surfaces of constant R. These 
are timelike within the locus f — and spacelike outside 
the locus. The values of z used to generate these curves are 
0.15,0.3,0.5 and 0.7. 




FIG. 6. The trajectories show surfaces of constant T. These 
are globally spacelike. The values of T used to generate these 
curves (bottom up) are tt/3, 7r/2,7r and 3tt/2. 



D. The fluid streamlines 



of course, globally timelike. From we have u > m/2 
for R > 2m. The scale factor is a(t) = C{1 - cos(T)) 1 / 3 
where C is a constant > 0. Writing e = r C/m we have 

R = (2 e (l-cos(r)) 1 / 3 + l) 2 

to 4e(l-cos(T)) 1 /3 ' ^ > 

and we note that R = 2m for T = To and T = 2ir — To 
where 

T = arccos I ^ 3 I . (33) 
Note that 1/2 < e < oo as explained in the Figure. 
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FIG. 7. Trajectories of constant r characterized by the con- 
stant e as explained in the text. We note that e — > 1/2 \/2 to 
the left and e — > oo to the right. 



VII. DISCUSSION 

Motivated by the AdS/CFT correspondence, we have 
examined the McVittie solution with a negative cosmo- 
logical constant A < 0. A detailed construction of the 
global structure has been given for the case of a back- 
ground of dust. We have found that the situation is very 
distinct from the cases A > pj . As in the pure Robert- 
son - Walker case, we find that A < ensures collapse to 
a Big Crunch, a feature which completely dominates the 
global structure. 
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The conformal representation of the fluid streamlines 
r = constant is shown in Figure[7] These trajectories are, 
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Appendix: x(T) 



x(T) = 



1 - 2.66T + 7.82T 3 

1.02 - 3.6T + 3.43T 2 - 1.9T 3 

0.557 - 0.772T - 0.395T 2 + 0.222T 3 

1.87 - 3.55T + 1.56T 2 - 0.235T 3 

-0.601 - 0.172T + 0.0217T 2 - 0.00230T 3 





0.177 
0.6 
1.42 
2.2 



< T < 

< T < 

< T < 

< T < 

< T < 



0.177 
0.6 
1.42 
2.2 

7T 



(A.l) 



